Time Reversal Symmetry breaking superconductivity 
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We study time reversal symmetry breaking superconductivity with Ak = A x 2_ y 2 (fc) + e lfl A Q (a = s 
or dxy) symmetries. It is shown that the behavior of such superconductors could be qualitatively 
different depending on the minor components (a) and its phase at lower temperatures. It is argued 
that such qualitatively different behaviors in thermal as well as in angular dependencies could be a 
source of consequences in transport and Josephson physics. Orthorhombicity is found to be a strong 
mechanism for mixed phase (in case of a = s). We show that due to electron correlation the order 
parameter is more like a pure d x 2_ y 2 symmetry near optimum doping. 



PACS Numbers: 74.72.-h, 74.25.Dw, 74.62.-c 

During last several years, a great deal of effort has 
been devoted to determine the symmetry of the super- 
conducting ordep. parameter (SCOP) of high T c cuprate 
superconductorso. While there is a strong concensus both 
from experiments and theory that the symmetry of the 
SCOP in the cuprates could be d x 2_ y 2 type; recent stud- 
ies provide increasing evidence that the pairing symme- 
try may be an admixture with a minor component, like 
s, d xy with predominant d x 2_ y 2 symmetry. The pairing 
symmetry provides clues to identify pairing mechanism 
which is essential for the development of the theory of 
high temperature superconductivity in cuprates, still re- 
mains a challenge after a decade of its discovery. 

Several different types of measurements which are sen- 
sitive to the phase of A(fc) indicate a significant mixing of 
s-wave component with d x i_ y 2 symmetry. For example, 
c-axis Josephson tunnelling studies on junctions consist- 
ing of a conventional s wave superconductor (Pb) and 
twinned or untwinned single crystals of YBCO indicate 
that the &C.OP of YBCO has mixed d±s (or d± is) 
symmetrjoB. Recently, a new class of c-axis Joseph- 
son tunneling experiments are reported by Kouznetsov 
et alp in which a conventional superconductor (Pb) was 
deposited across a single twin boundary of a YBCO sin- 
gle crystal. There measurements of critical current as a 
function of the magnitude and angle of a magnetic field 
applied in the plane of the junction provides a direct evi- 
dence for mixed d— and s-wave pairing in YBCO. A se- 
ries of high resolution measurements on thermal conduc- 
tivity (re) in the Bi^Sr^C aCu20% by Krishna et aZ,ushow 
that the re at low temperature becomes field-independent 
above a temperature dependent kink field Hk(T). This 
remarkable result indicates a phase transition separat- 
ing a low-field state where the thermal conductivity de- 
creases with increasing field and a high-field one where 
it is insensitive to applied magnetic field. The authors 
argue that this phase transition is not related to the vor- 
tex lattice because of the temperature dependence of the 
field Hk(T) (which is roughly proportional to T 2 ) as well 
as its magnitude. Instead, they suggest a field-induced 
electronic transition leading to a sudden vanishing of the 
quasi-particle contribution to the heat current. Possi- 
ble scenarios would be the induction of a minor e 10 A Q 



(a = d xy or s) component with d x 2_ y 2 symmetry with 
application of a weak field. A similar conclusion may 
also be obtained based on the results of angle resolve ph- 
toemission specroscopy (ARPES) experiment by J. Ma et 
al. in which a temperature dependent gap anisotropy in 
the oxygen- annealed Bi2Sr2CaCuOs+ x compound was 
found. The measured gaps along both high symmetry 
directions (T — M i.e, Cu-0 bond direction in real space 
and r — X i.e, diagonal to Cu-0 bond) are non-zero at 
lower temperatures and their ratio is strongly temper- 
ature dependent. This observed feature cannot be ex- 
plained within a simple d-wave scenario (but is consistent 
with d x 2_ y 2 -\-e l6 a scenario) and has been taken as a sig- 
nature of a two component order parameter, d x 2_ y 2 type 
close to T c and a mixture of both s and d othcrwisell. 
However, surprisingly enough, the experimental resultsEl 
was never been reproduced by any other equivalent ex- 
periments. In Josephson physics, it was foundEl that the 
current-phase relationship depends on the mutual ori- 
entation of the two coupled superconductors and their 
interface (the induced minor component then would be 
associated with any arbritrary phase, 0). This property 
is the basis of all the phase sensitive experiments probing 
SCOP. Therefore, motivated by strong experimental evi- 
dences of a mixed order parameter symmetry in high T c 
cuprates we present basic behaviors of superconductors 
with mixed order parameter symmetries such as d+e a, 
a = s, d xy , which is characterized by a local breakdown 
of time reversal symmetry. A clear evidence will emerge 
that the time reversal superconductivity at lower tem- 
peratures give rise to unusual behaviors depending on a 
and its phase. 

Assuming a well defined quasiparticle picture in the 
superconducting (SC) state, free energy of a supercon- 
ductor with arbritary pairing symmetry may be written 
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where Ek = \J (ife — fi) 2 + | A fc 2 are the energy eigen 
values of a Hamiltonian that describes superconductivity. 
It is generally believed that a single band Hubbard model 
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contains essential ingredients to describe Cu — O planes 
where transport processes takes place. Therefore, our 
model Hamiltonian reads as, 

"H = Hhub + Hp air (2) 

ija i i<7 

HpAin. = ^2v kk ,cl ] c t _ kl c_ k , l c k >> l (4) 

kk' 

where the notations have their usual meanings. 

Treatment of the Hubbard on-site correlation is a 
known problem. We use the Slave Boson theory of 
Kotliar and RuckensteirJain which the Hubbard model 
has solution for any filling and any value of repulsion. 
The si&ve - boson formulation of Kotliar and Ruckenstein 
(KR)B introduces four auxilliary boson fields correspond- 
ing to the occupancy of a site, namely empty (e^), doubly 
occupied (dj), singly occupied Pi- a ) with spin ±<r 

respectively. In terms of the slave boson field operators 
the single band Hubbard model takes the form 

Hhub = ]T Uj4Aa c 3<y z ]^ + U J2 d l d * - A* c ™ c ™ 

ij,(7 i ia 

+ ^2 A;(l - e\e t - d\di - ^p^Pia) 

i <J 

+ A ^(4t c *<t - d\di - p\ a p ia ) (5) 

where z la = (1 - d\di - p\ a Pta)~ 1,2 {e\pta +p\_ cr d i )(l - 

^i£i—p\- a Pi-o)~ 1 l' 1 and the Boson field operators are not 
independent of each other but are constrained by the re- 
quirements of completeness and local charge conservation 
at a site, hence the fourth and fifth terms are added in 
the Hamiltonian with the help of Lagrange multipliers (Ai 
and A-). The values of the boson field operators and La- 
grange multipliers are determined by minimizing the free 
energy in the saddle point approximation. This approx- 
imation diagonalizes the U term whereas renormalizes 
the hopping term as qt^ with q = z* z ; q usually being 
complicated function of u (= jj-), dopant concentration 
(5) and x(= e + d) carries the necessary information of 
electron correlation. In this approach solutions are ob- 
tained for thy paramagnatic state for all values of u and 
band fillingsEl that reproduces the correct Brinkman-Ricc 
result for metal-insulator transition (U c ) at half-filling. 

As a consequence the total Hamiltonian (2) takes the 
form, 

H = J2&k - + Y <fVk-'4T c -fci c -fc'i c fc'T 

ka kk' 

+ bosonic terms (6) 

where the band dispersion e k = 2»=i CiVi(k) are 
obtained from the angle rfesolved photoemission spec- 
troscopy (ARPES) results^ on Bi 2 Sr 2 CaCu0 8 . In de- 



riving Hamiltonian (6) fluctuations of the bose-field oper- 
ators are not considered (consequence of on-site correla- 
tion has shown up in (6) as e k — > qe k and V kk > — > V kk > = 
q 2 Vkk')- 

A mean field superconducting (SC) order parameter 
A* (As) = Ylk' ^ kk ' < c l:'-\ c -k'i > allows to solve the 
Hamiltonian (6) exactly to find the energy eigen values 
E k and hence minimization of the free energy (1) (i.e 
dF k ^ k i /dA k = 0) will result the superconducting (SC) 
gap equation as, 

^E^'^^^f 1 ), (7) 

where the pairing potential V kk > for a two component 
order parameter with a separable form and the corre- 
sponding gap functions are obtained as, 

Vkw = j2 v jfifi> & A fc = £y<*)A^ (8) 

3=0, j=0, 

For a mixed d x 2_ y 2 — d xy symmetry, Vo(i) = Vd aVS _ 2 

(Vd x „) and = cosk x ~ /3' cos k y (2 sin k x sin k y ) 

and the corresponding gap amplitudes are Aq^j = 
Ad 2 _ 2 (Ad xy ) (/?' < 1 corresponds to orthogonal or 
tetragonal symmetry). Similarly for a mixed d x 2_ y 2 — s 
symmetry, V\ = V s , fl= constant and Ai = A s . Substi- 
tuting equation (8) in (7) the gap equations for different 
components can be separated out as, 

k K 

The two coupled gap equations (9), a number conserving 
equation to fix chemical potential (/i) and a sixth order 
algebric equation in x(= e + d) (to get correct values of 
q for any value of u, S) are solved self-consistently for 
a fixed cut-off parameter Q c = 500 K around the Fermi 
level beyond which SC pairing does not exist. 

In figures l(a,b) phase diagrams of superconductors 
with A k — A x 2_ y 2(k) + e l6 a where a = A xy {k) (a) and 
=A S (b) are obtained for a fixed band filling p = 0.80 
and in a tetragonal square lattice. The phase diagrams 
comprises the amplitudes of pair condensates, evaluated 
at 5 K, in different channels as a function of the ratio 
7 = Vd 2 2 /V a . The minor components a appear at 
the same value of 7 irrespective of the phase 9 in the 
predominant d x i_ y i phase. But the rate of growth of the 
minor component is different depending on 9. For 9 = 
i.e for a real order parameter with (a) d x i_ y i + d xy or 
(b) d x 2_ y 2 + s wave symmetry, both the minor as well 
as the d x 2_ y 2 component increases with decrease in 7, 
thereby inducing each other. In case of d + s symmetry 
this phenomena is very pronunced. However, for any 
finite 9(^ 0) the d x 2_ y 2 component is suppressed with the 
enhancement in the isotropic s component at some values 
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of 7 depending on 9 resulting in a pure s-wave phase. 
This feature is absent in = A x 2_ y 2(k) + e 10 A xy (k) 
scenario. 
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cant role to the predominant d-wave. For example, when 
the minor component mixes with the <i-wave with a real 
coefficient {i.e., 9 — 0), the minor component has a very 
fast growth with lowering in temperature and it always 
induces the eZ-wave component by further stabilizing it. 
On the other hand, mixing of the minor component with 
a complex coefficient (i.e. 9 — 90°) suppresses the d- 
wave component when the minor component is s-wave 
and has no effect when the same is d xy . These behaviors 
are bound to have immense influence on thermodynamic 
quantities which are sensitive to the thermal gradient of 
the SC-gap. For instance, the specific heat will then have 
two jumps, at T c and at T" < T c where the minor com- 
ponent (a) vanishes and the nature of the second jump 
can be very different depending on the phase fl-providing 
signature of a true second transition below T^Ex 
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FIG. 1. Phase diagrams of superconductors with (a) Afc 
=A x 2_ y 2 (k) +e l0 A xy (k) and (b) A fc = A x2 _ y 2 (fc)+e l9 A a 
symmetries. Amplitudes of the d x 2_ y 2, d xy ~ and s— wave 
pairing condensates (at 5 K) are plotted as a function of 
the ratio of pairing strengths in the respective channels, 
7 = Vd 2 _ 2/Vdxy f° r u — 0-4. Note the qualitative differ- 
ence in (a) & (b); the onset of s-wave component influences 
the d a .2_j / 2-amplitude very strongly (suppresses or enhances 
depending on the phase 9) whereas the onset of d xy compo- 
nent have no influence on the d x 2_ y 2 -amplitude except for 
d x 2_ y 2 + d xy phase. 

For 7r/2 < 9 < 0, the minor component has two parts, 
real and complex (thereby giving rise to most probable 
mixing), that couples to the d x 2_ y 2 . A close look to figure 
1 therefore, would suggest that the mixing of a minor real 
s or d xy component will stabilize the d x 2_ y 2 symmetry, 
whereas inclusion of complex minor component s-wave 
tends to suppress d x 2_ y 2 gap. The d xy minor component 
with a complex phase has hardly any influence on the 
d x 2_ y 2 symmetry. These qualitatively different behaviors 
in different classes of mixed order parameter symmetries 
are also responsible for so with respect to thermal depen- 
dency. 

It is evident from figures 2(a,b) (for p — 0.8) that at 
lower temperatures the minor component plays a signifi- 
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Thermal behavior of superconductors with (a) 
y 2(k)+e l0 A xy (k) (7 = 0.195) and (b) A fc = 



FIG. 2 

A fc =A X 

A x 2_ y 2(k) +e' B A s (7 = 0.18) symmetries for various values 
of 9 (w = 0.4) . The minor (s or d xy ) component although 
have same T c their thermal growth inside the d- component 
is qualitatively different. Consequently, the predominant d 
component is influenced (e.g, further stabilised for 9 = and 
suppressed for 9 = 90) very differently. 

In figures 3(a,b) we present the angular correlation be- 
tween the amplitudes of different channels evaluated at 
5 K and p = 0.8. For any finite mixing of the minor 
component with d x 2 _ y 2 -wave thus produces states that 
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spontaneously breaks time reversal symmetry. Further 
interesting point to note is that with small orthorhom- 
bicity f}' = 0.98, (in the lattice as well), the s component 
gets enhanced largely (curves with solid circle symbol) 
whereas the change in d-component with orthorhombic- 
ity is rather dependent on 6. For example, the predom- 
inant (i-componcnt increases with orthorhombocity for 
< 6 < 0.227T or 1.787T < 6 < 2ir and decreases otherwise 
(both cases being small, cf. Fig 3(a)). The correspond- 
ing changes in Ad xy and associated change in 2 2 are 
negligibly small (cf. Fig. 3(b)) for such orthorhombicity. 
Therefore, orthorhombocity could be an intrinsic mech- 
anism for d&zs mixing (with an arbritary phase between 
them), d + s state being most favorable. 

In figures 4(a,b) we demonstrate the role of electron 
correlation on mixed pairing symmetry. It is shown that 
the minor component is always minimum at the optimum 
doping and is suppressed largely due to electron correla- 
tion explaining the novel features observed in ARPES 
experimentlij. This is based on the study of the ampli- 
tudes of the different components for the case of 9 = it/ 2 
(i.ed+is state) and for fixed T° = 45K at u = as a func- 
tion of density p (=1-8). The amplitudes are normalized 
with respect to the value of the d-component at optimum 
doping (the d-wave also gets suppressed with u but at a 
slower rate). Note, the electronic correlation suppresses 
the minor component very largely close to half-filling 





FIG. 3. Amplitudes of the d x 2_ y 2 and the minor gap (a) a 
or (b) d xy as a function of the phase of the minor component 
(for u = 0.4 & T" = 1AK). The curves with solid circles 
represent orthorhombic phase with (3' = 0.98 whereas the 
curves with open squares correspond to tetragonal lattice. 

whereas hardly have any effect at lower densities (the d- 
wave boundary shrinks a bit). Also, this behavior with 
density is similar for any other arbritrary value of 
Apparently, therefore, close to the optimium doping the 
order parameter will be more like a pure d x 2_ y 2 like. 
Furthermore, no mixing is found within a wide range 
of densities between d-wave with other probable minor 
components like s x 2 +y 2, s xy (e.g, these components have 
maxima at p ~ 0) 

Therefore, we have presented basic consequences of 
time reversal symmetry breaking superconductivity. It 
is demonstrated that mixing of a minor component at 
lower temperatures can give rise to unusual thermal be- 
haviors depending on the nature of the minor compo- 
nent as well as the phase associated with it. In the fully 
gapped phase at lower temperatures, as described by the 
Afc = A x 2_ y 2(k) + e l9 A a scenario, the number of quasi- 
particles will be exponentially small and hence will con- 
tribute negligibly small to thermal conductivity (and in 
any other transport properties). This scenario therefore 
can-|instify the experimental observation by Krishna et 
al.pB. This scenario also has strong potential to have 
natural explanation to the experimental observationu as 
discussed earlier, because mixing of a minor component 
(s or d xy ) with predominant d x 2_ y i wave would change 
sign across the twin boundary (hence the situation is de- 
scribed in figures 3) yielding a dramatically different an- 
gular dependence. 
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FIG. 4. Amplitudes of different components of (a) 
A fe = A x 2_ y 2(k) +iA s and (b) A fe = A x 2_ y 2(k) + iA xy 
superconductors at 5K. The minor (s or d xy ) component at 
u — has T™(0) = 45-?^. The electron correlation minimizes 
the mixing of the minor component at optimum doping. 
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In particular, this study provides a significant under- 
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ent phase. We also pointed out the subtle role of elec- 
tronic correlation with number density prohibiting mix- 
ing of the minor component at optimum doping. 
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